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Abstract 

We consider some cosmological aspects of nonlocal modified gravity 
with A term, where nonlocality is of the type RJ-(D)R. Using ansatz 
of the form OR = rR + s, we find a few a(t) nonsingular bounce cos- 
mological solutions for all three values of spatial curvature parameter 
k. We also discuss this modified gravity model from F(R) theory point 
of view. 



1 Introduction 

Modern theory of gravity is general theory of relativity, which was founded by 
Einstein at the end of 1915 and has been successfully confirmed for the Solar 
System. It is given by the Einstein equations of motion for gravitational field 
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5V : Rfj,v — \R9llv = 87rGT Mi/ , which can be derived from the Einstein-Hilbert 
action S = f y/—gRd 4 x + f y/—gC ma td 4 x, where g = det(g^ u ) and c = 1. 

Attempts to modify Einstein's theory of gravity started already at its very 
beginning and it was mainly motivated by investigation of possible mathemat- 
ical generalizations. During last decade there has been an intensive activity 
in gravity modification, motivated by discovery of accelerating expansion of 
the Universe, which has not yet generally accepted theoretical explanation. 
If Einstein's gravity is theory of gravity for the Universe as a whole then it 
has to be some new kind of matter with negative pressure, called dark energy, 
which is responsible for the accelerated Universe expansion. However, general 
relativity has not been verified at the cosmic scale (low curvature regime) and 
dark energy has not been confirmed in a laboratory. Hence, after discovery 
of the accelerated Universe there emerged a renewed interest in modification 
of Einstein's theory of gravity, which should be some kind of its generaliza- 
tion (for a recent review of various approaches, see pLj). However there is not 
a unique way how to modify the Einstein-Hilbert action. Among many ap- 
proaches there are two of them, which have attracted more interest than the 
others: 1) F(R) theories of gravity (for a review, see [21 13]) and 2) nonlocal 
gravities (see, e.g. [21 HI El El d El E] and references therein). 

In the case of F(R) gravity, the Ricci scalar R in the action is replaced by 
a function F(R). This has been extensively investigated for the various forms 
of function F(R). 

In the sequel we shall consider some cosmological aspects of a nonlocal 
gravity. Here, nonlocality means that Lagrangian contains an infinite number 
of space-time derivatives, i.e. derivatives up to an infinitive order in the form 
of d'Alembert operator □ which is argument of an analytic function. In string 
theory nonlocality emerges as a consequence of extendedness of strings. Since 
string theory contains gravity, as well as other kinds of interaction and matter, 
it is natural to expect nonlocality not only in the matter sector but also in 
geometrical sector of gravity. So the main motivation to consider an additional 
nonlocal term in the Einstein-Hilbert action comes from the string theory. On 
some developments in cosmology with nonlocality in the matter sector one can 
see, e.g., [101 EB H21 EES! E] and references therein. In the next sections we 
shall discuss a nonlocal modification of only geometry sector of gravity and its 
corresponding new cosmological solutions. 
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2 A Nonlocal Modification of Gravity 

In our consideration, nonlocal modification of gravity is a replacement of the 
Ricci curvature R in the Einstein- Hilbert action by a suitable function F(R, □), 
where □ = -^=d ll y/—gg ,JjU d v is d'Alembert-Beltrami operator. 

In this paper we consider nonlocal gravity model without matter, given by 
the action in the form 

oo 

where F(p) = ^^/n 1 ^™ is an analytic function of the d'Alembert-Beltrami 

n=0 

operator and C is a constant. Study of this model jT]) was proposed in [I] 
and some further developments are presented in [5J EJ [TJ [S]. This model is 
attractive because it is ghost free and has some nonsingular bounce solutions, 
which can solve the Big Bang cosmological singularity problem. 

By variation of the action ([!]) with respect to metric g^ v one obtains the 
corresponding equation of motion: 

c(2R^F(n)R - 2(V„V„ - giiU n)(F(n)R) - ^ v RF{u)R 

oo j, n—1 

+ E Y E fo" {g^d^Rd^-^R + O l RO n - l R) (2) 

71=1 1=0 

- 2d,a l Rd u a n - 1 - l R)^ = -^(G^ + Ag^). 
The trace and 00 component of ([2]) are: 

x /)— j 1 A 



Qn(F(n)R)+J2fnJ2( d v nlRd>lnn ^ lR + 2nlRnn ~ lR ) = ^r R - 



8itGC 2uGC 

n=l 1=0 

c(2R 00 F(n)R - 2(V V - g 00 n)(F(D)R) - l -g m RF(U)R 

oo „ n — 1 

+ E Y E (9oo (g aP d a riRdpU n - l - l R + U l RU n ~ l R) (4) 

71=1 1=0 

-2d n l Rd n n - l - l R)) = -±( Goo + Ag 00 ). 

/ 07TG 
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We use Friedmann-Lemaitre-Robertson- Walker (FLRW) metric ds 2 = —dt 2 + 
a2 {t) ( x^fc r 2 + T 2 d9 2 + r 2 sin 2 6d(p 2 ) and investigate all three possibilities for cur- 
vature parameter k (0, ±1). 



3 Ansatz and Solutions 

Investigation of equation (J2J) and finding its solutions is a very difficult task. 
In the case of the FLRW flat metric (k = 0) two nonsingular cosmological 

solutions for the scale factor are found: a(t) = ao cosh ^<y/|^, see [HE], and 

a(t) = aoe^v^f* 2 , see [7]. 

To get some new solutions we also use ansatz of the form 

DR = rR + s, (5) 

proposed in [5j, where r and s are real parameters that will be fixed later. The 
first two consequences of this ansatz are 



n n R = r n (R+-),n>l, JF{U)R = F(r)R+ -(JF( r ) - /„). (6) 

Now we can search for a solution of the scale factor a(t) in the form of a 
linear combination of e xt and e~ xt , i.e. 

a(t) = a (cre A * + re~ xt ), < a , A, a, t e R. (7) 

Then the corresponding expressions for the Hubble parameter H(t) = -, scalar 
curvature R(t) = -^(a'd + a 2 + k) and OR are: 



H(t) = 


\{ae xt - Te~ xt ) 
ae xt + re~ xt ' 


R(t) = 


6 (2a 2 X 2 (a 2 e ux + r 2 ) + ke 2tx ) 


a 2 (ae 2tx + rf 


OR = 


12X 2 e 2tx {Aa 2 \ 2 aT - k) 


a 2 {ae 2tx + r) 2 


cLS 





(8) 



OR = 2X 2 R - 24A 4 , r = 2A 2 , s = -24A 4 . (9) 



4 



Substituting parameters r and s from (jSJ) into (J6]) we obtain 

□*\R = (2A 2 )%R-12A 2 ), n>l, 
F(n)R = F(2\ 2 )R- 12A 2 (^(2A 2 ) -/„). 1 j 

Using this in (j3J) and (j4j) we obtain 

36A 2 J^(2A 2 )(i? - 12A 2 ) + T (2A 2 ) (4A 2 (i? - 12A 2 ) 2 - i? 2 
-24A 2 / (^-12A 2 ) = |^, (11) 

(2i? 00 + X -R) [7{2\ 2 )R - 12A 2 (^(2A 2 ) - / )) - ^'(2A 2 ) [r 2 + 2A 2 (i? - 12A 2 ) 2 

- 6A 2 (^(2A 2 ) - f )(R - 12A 2 ) + QHF(2\ 2 )R = -^L_(G 00 - A). (12) 

Substituting a(t) from (j7j) into equations (TTTT) and (TT2"|) one obtains respec- 
tively the following two equations as polynomials in e 2Xt : 

( 3A 2 _ A) + 3a 2 r 4 Q ie 2A < + 6a 2 ar 3 Q 2 e 4At - 2arQ 3 e 6Xt 
+ 6a 2 a 3 rQ 2 e 8A ' + 3a 2 a 4 Q ie 10A ' + ^ (3A 2 - A) e 12A * = 0, (13) 



8ttG 



£ (3A 2 - A) + 3r 4 a 2 i? ie 2A * + 3r 2 R 2 e 4Xt + 2arR 3 e e 



3a 2 R 2 e 8Xt + 3a i a 2 R l e wxt + (3A 2 - A) e 12At = 0, (14) 

87rG 



where 



Qi = 36CA 2 i^(2A 2 ) + a 2 (-96C/ A 4 + - -^)ar 

7T(_t Z7TG 

+ 24C/ A;A 2 + ^, (15) 

7A 2 5 A 

Q 2 = 72CX 2 KT(2X 2 ) + a 2 (-192C/ A 4 + ^ - ^)<rr 

+ 48C/ofcA 2 + i ^, (16) 
Q 3 = -324Ca 2 A 2 ariG^(2A 2 ) + 144CA 2 K 2 J r/ (2A 2 ) 

- a 2 M216C/ A 2 + ^)ar + a 4 (864C/ A 4 - ^ + (17) 
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3A 2 -A „ 

Ri = Qi- A%G oral, (18) 
R 2 = _6C (k - 12a 2 X 2 aT) KF{2\ 2 ) - 36CA 2 K 2 J r, (2A 2 ) 

+ (1927tGC/ A 2 + 1) ar - (3072ttGC/ A 4 + A 2 + 5A) a 2 r 2 , 

LllKj o7rGr 

(19) 

#3 = -18C (fc - 6ajjAVr) fCF(2A 2 ) + 36CA 2 K 2 J r '(2A 2 ) 

+ ^ (192^GC/ A 2 + 1) ar - (34567rGC/ A 4 + 3A 2 + 5A) a 2 r 2 , 

(20) 

and X = 4aoA 2 crr — k. 

Equations ( !T3|) and (TT4"]) are satisfied when A = ±yf > as well as Qi = Q2 = 
Q3 = and i?i = i?2 = R3 = 0. Note that this approach to find conditions 
under which solution exists differs with respect to approach used in [5] and [TJ. 

The corresponding solutions can split into the following three cases. 

Case 1. 

*(2A")=0, *<(2tf)=0, /o^-^. (2D 

Case 2. 

3A; = 4a 2 Aar. (22) 

Case 3. 

^ (2A2) = %^CA + ^ /o ' r ( 2A2 )= ' fc = -4aSAar. (23) 

In the first case we have family of solutions for arbitrary a, r and ao 

a(t) = a (ae xt + re~ xt ) 

with function J 7 satisfying conditions given in ( !2~TT) and arbitrary A; = 0,±1. 
Thus it includes also solution a(t) = a cosh C\J~^t)i which was found in [4] 
with addition of some radiation [5] in the action. 

The second case yields a family of solutions for arbitrary a 7^ and ao 

a(t) = a ( ae xt + e~ xt 
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which are valid for arbitrary analytic function F . 
The third case yields another family of solutions 



and function F has to satisfy conditions given in ( |23i) . 

Note that if k = then equation ( 122]) and third equation in ( 1231) coincide, 
cr or r has to be zero, and we have 2 solutions 

a x {t) = a e xt , a 2 {t) = a e~ xt , 

where a is absorbed into ao- These are the de Sitter solutions, see also [8]. 



4 F(R) Theory 

Action ([T]) with ansatz (jSJ) gives an theory and we get another way to 

analyze the above solutions . The corresponding action for F(R) is 

S' = I ^L-F(R)d 4 x, F(R) = aR 2 +/3R- 2A, (24) 
J lbnG 

a = 8nGCF{2\ 2 ), f3 = 1 - 96ttGC\ 2 (F(2\ 2 ) - f ). (25) 
By variation of the action (|24p we get the following equation [2]: 

F'(R)R lll/ - \F{R)g^ - [V M V, - g^n]F'(R) = 8nGT^. (26) 
The generalized Friedmann equations are 

/d N 2 k 8ttG , 

y + ^=w(R) (p+ ^ 

a = -3FW (P + 3p + 7J + 3p) - 

From the equations of motion we derive expressions for effective density p 
and pressure p 

P = ^ (~F(R) - ^F'(R) + 3±RF"(R)) , 

87tGV2 2 a / ^7) 

P=^a &(R) ~ §F'(R) - (3^ + R)F"(R) - R 2 F"\R)) . 



In the case of action S' we have 



1 (a. 9 a ■ 

P =8^(f fl! + A + 2a(3 ? + ^)' (28) 
w = - = -l V a 



P f R 2 + A - 6a§ R 

It seems to be natural that any F(R) theory should satisfy [2] 

F'(R) > 0, -1.1 < w < -0.98. (29) 

In the Case 1, we have J-"(2A 2 ) = and consequently a = and therefore 
F'(R) = p = §, w = -1. 

In the Case 2, we have a solution of the form a(t) = ao(ae xt + re~ A *), where 
the constants a , cr, r, A are subject to the constraints 

A 2 = -A, 3k = Aether. (30) 

From this set of parameters we see that 

R(t)=4A, w = -l. (31) 
In the Case 3, we obtain 

a= ik + ^ir^ fo > / 3 = ^( 1 + 16 ^ A /o), (32) 



and 



where 



w = —1 + h(k, u, v), (33) 



512ku {k 2 + I2ku + 144n 2 ) (v + 1) 

{ > U ' V) - (jfe4 + 20736tt 4 ) Vi + (48fc% + 6912&U 3 ) V 2 + 288£rW 3 ' 

Vl = 2u + 5, 1/ 2 = 14v + ll, V 3 = 70^ + 79, (35) 

m = -agAaWr*, « = 647rGC/ A. (36) 

Obviously, in the case k = we have w = —1. Note that h{— l,u, v) = 
h(l, —u,v). Hence it is sufficient to analyze only h(l,u,v) and it will be pre- 
sented elsewhere with some other properties of this F(R) theory. 



8 



5 Discussion and Concluding Remarks 



In this paper we have considered a nonlocal gravity model with cosmological 
constant A and without matter. We found three types of nonsingular bouncing 
solutions for cosmological scale factor in the form a(t) = ao(ae xt + re~ xt ). 
Solutions exist for all three values of spatial curvature constant k = 0, ±1. All 
these solutions depend on cosmological constant A, which is here an arbitrary 
positive parameter. 

Note that trace equation (ITTj) can be rewritten in the form 



A 1 R + A 2 (4A i? - R ) + A. 



0, 



where 



A, 



1 



12A 2 (3J r (2A 2 ) - 2/ ) - 96A 4 J r/ (2A 2 



8nGC 
A 2 = r(2X 2 ), 

A s = ^4t^ - 144A 4 (3j^(2A 2 ) - 2f ) + 57QX 6 J r '(2X 2 



From A\ = A2 



2ttGC 

A3 = one obtains the following system of equations: 

1 



12A 2 (3^(2A 2 )-2/ ) 
144A 4 (3J r (2A 2 ) - 2/o 



8ttGC" 
A 

= 2-kGC 



J z '{2\ 1 ) = 0, 



Finally, we get 
A 2 



A 
3"' 



(37) 

(38) 
(39) 
(40) 

(41) 
(42) 

(43) 



and this is related to the corresponding conditions in [5] and to the first two 
equations of our Case 3. However in Case 3 there is additional condition 
k = — 4a 2 1 A<7T, which does not permit solution of hyperbolic cosine type when 
k — 0. In [H [5] this problem was solved adding some radiation in 00 equation 
of motion. 

Note that in the above solutions for the scale factor one can write more 
general expression by replacement t — > t — to, i.e. 



a(t) 



a (ae x{t - to) +re- x{t - t " ) ). 



(44) 
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When a > and r > 0, and to — 2vln(— ) then ( I44p can be rewritten as 

a(t) = 2a ^/aT cosh(Xt) . (45) 

If t < < a and t = ^-ln(— -) then the scale factor (j4"4"|) becomes 

a(t) = 2a 0V / -(7T sinh(At). (46) 

Note that one can construct a solution a(t) for a negative value of A. In 
this article we used linear ansatz fl5]), but some other ansatze [15] can be also 
useful in search of new cosmological solutions. Here we considered a nonlocal 
gravity model without matter, however modern cosmology will probably need 
both - modified gravity and an exotic matter (see, e.g. zeta strings [16] and 
p-adic matter [T7]). 
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